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Abstract 

Isothermic parameterizations axe synonyms of isothermal curvature 
line parameterizations, for surfaces immersed in Euclidean spaces. We 
provide a method of constructing isothermic coordinate charts on sur- 
faces which admit tlicm, starting from an arbitrary chart. One of the 
primary applications of this work consists of numerical algorithms for 
surface visualization. 

MSG 2010: 53A10 

Keywords: isothermal coordinates, isothermic coordinates, isothermic 
surfaces. 



1 Constructing Isothermic Coordinates (Curva- 
ture Line Isothermal Charts) 

Isothermic parameterizations are synonyms of isothermal curvature line pa- 
rameterizations [21 01 |5], for surfaces immersed in Euclidean spaces. By 
definition, isothermic surfaces represent immersed surfaces which admit at 
least one isothermic parameterization. The initial motivation of this project 
came forward from solving Bonnet problems analytically and numerically, by 
using Cartan's theory of structure equations, and numerical analysis. The 
authors also studied dual surfaces, in the sense of Christoffel transforms. 
Such transforms exist exclusively for isothermic coordinate charts. 

In a previous work [T], we solved Bonnet problems in terms of Car- 
tan theory in conformal coordinates in the most general context of regular 
surfaces. That is, for any compatible first and second fundamental form, 
corresponding families of solutions were obtained (up to roto-translations) . 

Due to the high impact of isothermic coordinates in various areas, it 
became important to study necessary and sufficient conditions for the exis- 
tence of isothermic coordinates on surfaces; and whenever such coordinates 
exist, a practical method of constructing them became imperative. 

We provide a method of constructing isothermic coordinate charts on 
surfaces which admit them, starting from an arbitrary chart (parameteri- 
zation). One of the primary applications of this work consists of numerical 
algorithms for surface visualization. Many of these results can be generalized 
to immersions of 2-dimensional manifolds in n-dimensional space forms, but 
our main goal was a construction algorithm with numerical implementation 
and visualization. 

2 Fundamental Concepts of Surface Theory 

This work exclusively involves immersed (respectively, embedded) surfaces 
in the Euclidean 3-dimensional space. 

Definition 2.1. A local surface is a differentiahle mapping from an open, 
simply connected subset of the Euclidean plane, to M'^. Various references 
use the name of patch of a surface, or local patch instead of local surface. 
We will call such a local surface a surface immersion, if the Jacohian of 
the differential mapping has rank two at every point. This is equivalent to 
the differential map being 1-1. Remark that such a surface is generally not 
embedded. If the local surface (immersion) represents an injective mapping 
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in itself, then we will call such a mapping an embedding, or an injective 
patch. 

In a more global sense, a connected set M in is said to be a regular 
2-dimensional surface if for an arbitrary point p & M there exist an open 
ball Up in with center at p and a homeomorphism ip which maps M DUp 
onto an open disk on some plane in the space R"'. Such a homeomorphism 
is frequently called a chart, or a system of local coordinates on a surface. 

It is worth noting that, in the previous definition, the word disk can be 
replaced by an arbitrary open set of a plane, diffeomorphic to a disk. 

Let us consider an immersion r = r{u,v) = {X{u,v),Y[u,v), Z{u,v)) 
(of an open, simply connected set) in M^. Note that the property of this 
map r = r{u,v) being an immersion is equivalent to the property that the 
vectors ^ and ^ are linearly independent at every point. These span a 
tangent plane at each point, that is, generate a tangent bundle. We will 
denote N = i''"^^", and call it the unit normal vector field, or Gauss map; 
it represents a map N : D ^ S'^. 

Definition 2.2. We will call the following quadratic form 

ds'^{u, v) =< ru, ru > du^ + 2 < r„, r^, > dudv+ < r^, > dv"^ 
a naturally induced metric, or first fundamental form. 

Recall that the notion of induced metric came from a bilinear map that 
acts as follows: I{vp, Wp) =< Vp, Wp >, where the vectors Vp and Wp represent 
tangent vectors to the surface at the point p. Also recall that, locally, every 
smooth surface can be represented as a graph (Monge chart). In general, 
we will not work with Monge charts, but with generic immersions. 

If we define E = Ir^P, F =< ru,rv >, G = |r„p, then we can identify 
the naturally induced metric (2-form) 

ds^ = Edu^ + 2Fdudv + Gdv^ , 

with a symmetric and positive definite real-valued matrix. 

As it is very well known, the notations above are classical, and were first 
used by Gauss. 

Definition 2.3. For a local immersion in M"^, with image M = r{D), let N 
represent the Gauss map. For a tangent vector v = Vp to M at p, we put 
Sp(v) = -dNpiv) = -D^N{p). 

Then Sp, viewed as a linear map from the tangent plane to itself, defines 
the shape operator S at every point. 
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We define the normal curvature as 

< Vp, Vp > 

Clearly, if Vp is unitary, then the denominator is equal to 1. 

Note that, if the actual immersion a surface is unknown and only the first 
fundamental form is given, then we do not possess much information about 
the geometric properties of the surface. By a classical result of Bonnet, the 

geometry of surfaces exclusively depends on two quadratic differential forms. 
One of them is the metric, and the other is the second fundamental form. 

Definition 2.4. The second fundamental form of a regular surface M in 
represents a symmetric and bilinear form on the tangent plane to M at 
the point p given by: 

II{vp,Wp) =< S{vp),Wp > for any arbitrary tangent vectors on the tan- 
gent plane at p. 

Another classical notation of the second fundamental form that is due 
to Gauss is: 

// := Idu^ + 2mdudv + ndv^ , 

where the coefficients I, m, n are given by: I =< rum N >, m =< ruv, N >, 
n =< rvv, N >. 

Unlike the first fundamental form, the second fundamental form is not 
necessarily positive definite. However, it is a symmetric bilinear form. Since 
all directions in a tangent plane TpM form a compact set homeomorphic 
to a circle, k has at least one minimum and one maximum, i.e., at least 
two extremal values. The maximum normal curvature ni and the minimum 
normal curvature K2 are called principal curvatures. They represent the 
eigenvalues of the shape operator, and the eigenvectors corresponding to 
them arc called principal directions. 

Definition 2.5. The average and the product of the two principal curvatures 
are called mean curvature and Gaussian curvature, respectively. 

In the classical literature, surfaces whose mean curvature (respectively, 
Gaussian curvature) is a constant have been studied extensively. Constant 
mean curvature (CMC) surfaces represent a special subclass of Bonnet sur- 
faces, which in their turn represent an important subclass of isothermic 
surfaces. 
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Definition 2.6. // the induced metric is a (nonconstant) multiple of the 
flat metric, then the immersion is said to be conformal, or isothermal, or in 
conformal coordinates, i.e. \\ru\\ = and < ru,r^ >= at every point. 

A surface is said to be (parameterized) in curvature line coordinates, if 
both first and second fundamental forms are diagonal. 

A surface in M'^ is said to be parameterized in isothermic coordinates if 
the given parameterization is conformal and in curvature line coordinates at 
the same time, i.e. it is given as an immersion f : D <Z B? — > E? such 
that 

< fujv >= 0, 

and 

\\fJ^ = \\fJ^ = E{u,v), 
for a smooth function over D, and < fuv, N >= 0. 

Remark that obtaining a parameterization that is both isothermal and in 
curvature hne coordinates at the same time is not possible on most regular 
surfaces. 

Also note that at umbilic points (i.e., points where the normal curvature 
is the same in every direction), the actual choice of an orthogonal system 
of eigenvectors becomes problematic. Therefore, we usually view all the 
umbilics as singularities. Clearly, this point of view makes the sphere a very 
special example of isothermic surface! 

Assume that we are dealing with a regular surface in which admits 
isothermic coordinates (away from a discrete set of singularities). We will 
call such a surface an isothermic surface. 

Isothermic surfaces were extensively studied in [6] , where the main inter- 
est was to study Bonnet surfaces as a particular case of isothermic surfaces. 

It is important to note that |6] mentioned the following important char- 
acteristic property of isothermic surfaces: 

M is an isothermic surface in if and only if locally there exists a 
conformal parameter z = u + iv such that the Hopf differential coefficient, 
Q = Q{z,z) is a real-valued function. 

This characteristic property (in terms of Q being real) is obvious from 
the definition of an isothermic surface, together with the formula of Q in 
terms oil, m and n. 

A more comprehensive characterization of isothermic surfaces was pro- 
vided by Bobenko in [3]: 
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Lemma 2.1. Let f : D — > be a conformal immersion of an umbilic-free 
surface. The surface is isothermic if and only if there exists a holomorphic 
non-vanishing differential f{z)dz'^ on D and a function q from D to 
such that the Hopf differential is of the form: Q{z,z) = f(z)q{z,z). It is 
easy to see that, in this case, w = J ^/ f{z)dz is an isothermic coordinate. 

We arrived at the following question, which is important to answer in 
this work, from both a theoretical and a constructive (applicative) point of 
view: 

Question: "Starting from a arbitrarily given immersion {x, y) i— )■ f{x, y) 
of an isothermic surface in the Euclidean 3-space, is there a simple method 
of obtaining an isothermic parameterization i->- f{P,j) corresponding 

to it?" 

We implemented a key idea for changing the old coordinates to the 
isothermic ones, which can be summarized in two steps: 

a) , applying the Gram-Schmidt orthogonalization; 
and 

b) . multiplying the velocity vectors by a 'non-constant scaling func- 
tion' K = K{x,y), namely a class non-constant function specific to this 
chapter (not to be confused with the Gaussian curvature). 

Corresponding to this method, we obtained a condition that this func- 
tion, K, must satisfy, and we constructed isothermic coordinates from an 
arbitrary immersion, under the above-stated condition. 

Solving a Bonnet problem gives rise to a generic immersion, which is 
rarely isothermic; constructing an isothermic parameterization for the same 
physical surface, whenever possible, presents a great deal of simplification. 

An isothermic parameterization creates a special mesh that is desirable 
in some applications, such as: 

• applying geometric methods to mesh discretization which are both 
homogeneous and orthogonal; 

• in structural mechanics, deriving the equation of shells along the prin- 
cipal directions; 

• in architecture/industrial design, working with orthogonal meshes is 
often preferrable, in order to minimize the production cost; 

• finding dual surfaces (Christoffel transforms). 
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Theorem 2.1. Let f{x,y) = {f'^{x,y), f^{x,y), f{x,y)) be an arbitrary 
immersion, with velocities fx and fy, of a regular surface M = f{D), where 
D represents an open and simply connected domain. Let us define the vector 
distributions 

/, := K(.. y) c„.(a) + ^f^^ s.n(a)) , (2.1) 

/. := K(., y) (- -| sMa) + ^|^^ cos(a)) , (2.2) 

where the function K{x, y) is positive and smooth, and the angle a satis- 
fies the following equation (in terms of coefficients of the first and second 
fundamental forms corresponding to f): 

Further assume that there exists a change of coordinate map {x, y) 
7) such that the given vectors fi and /2 respectively represent the partial 
velocities of f . 

In these assumptions, the new coordinates 7) are isothermic 
if and only if the corresponding smooth function K, with respect to 
the two different charts, will satisfy the following equations: 



2 / 'E^Fx "2 J7]Fx ~I~ Gx , , Fy J^Fx ~\~ 2-^2; _ 



Kj = K 



2(G-f )v^ 



cos (a) — 



2fJg--4 



sin(a) 



_ F 

y 



sinfa) + , 



El 

E 



cos(a) , 



(2.4) 



E — 2—F -4-0 



2{G 



F 



sin(a) — 



y E ^ I X 



cos a 



2fJG-^ 



cos(q) 



_ F 

"y E^^ 



G 



i?2 



sin(a) 



(2.5) 
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Remarks about the previous theorem: 

A) . Note that the geometric meaning of a{x,y) is of a rotation angle 
that changes at every point (of a pre-orthonormahzed frame after Gramm- 
Schmitt); and the geometric meaning of K{x, y) is of a non-constant dilating 
(scaling) function of the unit tangent vectors, up to a pair of partial veloc- 
ity vectors in isothermic coordinates at each point on the surface (both of 
magnitude K{x,y)). 

B) . Consider the fractional expression of tan(2a). If the initial param- 
eterization would be an isothermic one, one can easily check that both the 
numerator and the denominator would become zero at umbilic points. For 
umbilics, the expression would then be undetermined, and not undefined. 
In the following statement we assume that we are always away from umbilic 
points. However, even for the case of isolated umbilic points, the discrete 
algorithm allows a certain flexibility in choosing the size steps in such a way 
that the umbilic points can be avoided. 

C) . The formulas above provide a simple method to reparameterize to 
isothermic coordinates. We are not claiming that this is the only method. 



D). In this theorem, equations (2.1)-(2.2) and (2.4)-(^2.4| assure the com- 



patibility condition (/i)/? = (72)7, which will become: (/'y)/3 = (7/3)7. 

Proof. Note that our long expressions in the proof will not be simplified to 
the maximum extent, due to our willingness to show all the computations 
exactly as they were successively performed. 

By the definition of isothermic coordinates, the vectors /i = and 
/2 = represent the partial velocity vector fields of an isothermic parametriza- 
tion (/3, 7) if and only if they satisfy the following conditions simultaneously: 

1- 11/711 = 

2. < fj, >= 0, 

3. <f^p,N>=^. 

4. U,)p = (//3)7- 

1. It is easy to rewrite the first condition in an equivalent way. Specifi- 
cally, we obtain 

ll-fl|2 ^ r r ^ r^2|ll/a:lP 2 / \ , f ^ ^ f V ^ ^ IsWf I \ ■ I \\ 

WfiW =< f-i! fj >= COS (a) + 2 cos(Q)sm(a) 



ll/.ll^ + i^l/.'" 



E 
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and 



m\' =< U,U >= te£sin^(a) -2lAlA^^^cos(a)sin(a) 



G- 



E 



Note that the ||/-y|| = ||/y3|| = K is an immediate consequence of the 
definition of the initial vector field distribution fi,f2- 

2. The second condition is also immediately implied by the two vector 
fields provided 

, J. _f . I I I I I ^ / \ ■ / \ , f^l fy ~E\\f^\\ 2 / N 1 

< 77 ) //9 >= K - ' ' cos(q) sin(Q) H , cos (a) 

, 7V-2 / fx, fy > — -gWfxW . 2 / X \ 

+K sm (q) 



G- 



El 

E 



+7^ I ^ sm(Q) cos(a) 



E 



3. We will show that even the third condition is also implied by the vector 
construction in the theorem hypothesis. 

In order to reformulate the third condition, let Di{x,y) = and 
D2{x,y) = ^j=J=, for simplification. 

Therefore, fj and fp will become 

:= K {Di cos(a) + D2 sm(a)) , (2.6) 
:= K {-Di sin(a) + D2 cos(a)) . (2.7) 
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For the last identity, consider 



— D 

< if-y)0,N > = < K^iDi cos(a) + D2 sin(a)) + {- sm{a) cos(a) 



9 n — 



cos [aj H sm(a) cos(Qf) 



p2 



G- 



p-2 



■ — — °^ COS (a) H sm (q) — — — ° sin(a) cos(a) 



- sm(Q) cos(q) sm(Q) cos(q) + 



G 



G 



p2 



cos^{a) 



: sm(a) cos(a) — — 



F Dia^ 



G- 



G-^ 



: COS (a)), TV > 



(2.8) 



Since < fx, N >= and < fy, N >= 0, then aU terms with < I^i, iV > 
and < D2 , > will be z ero; replacing the partial derivatives of Di 
and D2 into the equation (2.8) we obtain 



(cos^(a) - sin^(a)) < fxx,N > 



+ 



E{G-§ 



i(a)cos(a) < fxx;N > 



+ 



+ 



(cos (a) — sin (a)) — 2 ^ ^ sin(Q) cos(a) | < f^y, N > 



G-^VE 



G 



G- 



p2 



cos(a) sin(a) < fyy,N > 



F 
' E 



G-^^/E JG- 



-1 + 



p2 ^ 



-m I cos(2a) 



+ 



G 



EG 



1-2- 



G- 



p2 



m + 



sin(2a) 



G-^E 



(2.9) 



Now, if we rewrite the equation (2.3) as 
/ 

sm(2a) = -2 



2^-G 



1-2- 



f2 



m + 



f2 



cos(2a), 



10 



and plug it into the equation (2.9), then we obtain 

1 



< {f-,)p,N>-- 



F 
' E 



-1 + 



m cos(2a) 



- G 



ElG 



p2 



1-2 



1 + 



m + 



f2 



cos(2q) = 0. 



-1-2 



\ 7 ir'i \ I' ^ ^ H ^PT I 



4. Last but not least, we will analyze the fourth condition in its equivalent 
forms, and show it is actually equivalent to the PDE system satisfied 



by the function K. Again, consider the equations (2.6)-(2.7) 



(/^),3 = Kfi (Di cos(a) + D2 sin(a)) + K— (Di cos(a) + D2 sm{a)) 



, , T?-(-Di cos(q) + D2 sinfa)) 
Kp{Di cos(a) + D2 sm(a)) + K{~-S2^ ^ ' ^ 



^(Di cos(q) + -D2sin(Q)) - f cos(a) + ^Da sin(a)) 



G 



F2 



KpiDi cos(a) + D2 sin(a)) + sm(Q) cos(a) - sin^(a) 



sin(Qf) 
cos(q)) 



'E 



+ 



dy 



F2 



F i^Di 2 
" sin(a) cos(a) — — ° cos (a) 

„ p2 E _f2 

- "B V *^ ^ 



+ sm (q) '^'^ 



cos(a) 



G 



-D2Q3: 
Vi5 



i(a) cos(q) 



G- 



F2 



: sin(a) cos(a) + 



D20ly 



G 



p2 



: COS (q) 



, F Dia^ F D2a^ 2, 

+ — — , sin(Q) cos(q) — — — cos (a)), 



E I^_P2 
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d 

{ff})j = Kpi-Di sin(a) + Da cos(a)) + K — {-Di sin(a) + D2 cos(a)) 

, , sin(a) + 7)2 cos(q)) , , 
= Kp{-Di sin(a) + D2 cos(a)) + ^= ^ cos(a) 

|-(-Disin(a) + D2Cos(a))- f |:(-Disin(Q) + D2Cos(a)) 
+^ / sin(a)) 

— D —D 
= K-,{-Di sin(a) + D2 cos(q)) + Js:^(-^3lJ: sin(a) cos(a) + cos^(a) 

. 2, , , ■ , , , , , F i^Di . 2, , 

: sm (a) H sin(a) cos(a) + — ^ sm (a) 



-D1Q3; gjjj2 / -> _ r — sin(a) cos(a) — ""'^'^ sin(a) cos(a) 
fE ^ ' E _ Iq_ 



DiOy . , , , , -DaOj, . 2. X 
sm(aj cos(Qj . sm (a) 



+ F / = sin(Q) cos(q) + - ^== Sin (a)). 

The compatibility condition is rewritten as 

= (/7)/3 - (/,3)7 = ^f/3(-Di cos(a) + D2 sin(a)) - K^{-Dx sin(a) + 1)2 cos(a)) 

^,D2 F jl-Di L>ia, D2as; F -D2Q. 



Kp{Di cos(a) + D2 sin(a)) — K^{—Di sm(Q) + D2 cos(q)) 



^pr2^ fyx- %fxx- -§^{i)fx Uv Ef^^dxiyG fi) 



+ 



Diax , D2ay F D2ax 



^ s/g^ ''sIg^ 

After some simplifications, we get a rephrasing on the compatibility 
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condition (/i)^ = (7^)7 in the equivalent form 

= {fi)(i ~ {ffi)i = {Kfi cos(q) + sin(a)) Di 



+ I Kp sin(a) - cos(a) + \ + ^ ' 



I?2 



Kfi cos(a) + K^i sin(a) + K 



-Ey — -jjiEx + 2Fx ^ ax 



2Ex/G-^ 



'E 



+ \Kp sin(a) - K., cos(q) I ^^^^ 2 ^ + G, 



2(G-f 



+ 



F2 



G 



^2 



(2.10) 



Since the compatibiHty condition above has to be satisfied at every 
point of the domain and since Di and D2 are hnearly independent, 
the condition 4) is equivalent to the following PDE system 



Kp cos(q:) + K-y sin(Q) — K I — 



2 1 Fy Ex ^Fx 



2eJg^^ 



E 



(2.11) 



Kp sin(a) — K^ cos(a) = K 



']^Ex — 2-^Fx + Gx 



2(G 



/E 



n f2 



(2.12) 



Remark however, that the previous system (2.11), (2.12) in Kp and 
is compatible and has a unique solution, namely the PDE system 



(2.4) satisfied by K in the statement of the theorem. 



□ 



On the construction algorithm of the isothermic parameterization (/3, 7) 
starting from the initial, arbitrary parameterization {x,y), it becomes nat- 
ural to write down the actual formulas which describe this transformation 
of coordinates (change of chart). The Jacobian of this change of chart can 
be easily computed, and is presented in the following Proposition. 

Corollary 2.1. The Jacobian of the change of chart 7) 1— )■ {x,y) has the 
following entries: 



^ 'cos(a) F sin(Q!) 



xp = K 



— sin(Q;) F cos(a) 



(2.13) 
(2.14) 
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and 




(2.15) 



(2.16) 



Moreover, each of the PDE systems {x,y) listed above, namely (2.13)-(2.14) 
and ( |2.15 )-(2.16), has a unique solution, given a choice of initial values 
{xo,yo). 



Proof. It is straightforward to verify that the compatibihty conditions, x-yj3 = 
xp^ and y^p = yp^, are verified, because the condition f^j^ = fjj^ is already 
fulfilled. Given the initial value 3;(0,0) and y(0, 0) by results from [7J, the 
PDE system above has an unique solution represented by the chart (x, y) 
which corresponds to (/3,7). □ 

In the previous theorem, we proved that the existence of a certain scaling 
function K{x, y) of arbitrarily given coordinates is necessary and sufficient 
in producing some new, isothermic coordinates. K represents the solution 



of a system of differential equations, (2.4). Moreover, we can provide a 



condition for such a scaling function K to exist and be unique. 



Theorem 2.2. As in the statement of the previous theorem, let E, F, G and 

l,m,n represent the coefficients of first and second fundamental forms of the 



arbitrary immersion f{x,y), and let a be given by the same condition (2.3) 
as before. Assume the change of chart (x, y) i— )■ (/3, 7) is introduced in the 
same way as in the previous theorem. Let us now consider the PDE system 



as defined by (2.4). 

The following condition 



/a. 

dy 



2{G- 



E dx 



- E —0 — E 



n F^ 



2(G- 



'-)s/E 



J5_ 
dx 



+ 



+ 



F± 



2(G- 



9 ( ax \ 
dx \ ^) 



'E 



+ 




dy 



+ 



4^ 



) + (T7fcK) 



E 



= 0, 
(2.17) 
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assures the existence of a function K{x,y) that satisfies the system of dif- 
ferential equations (2.4). Moreover, if the solution K exists, then it will he 
unique, given some appropriate initial value condition. 



Important Remark. The uniqueness of K that satisfies the PDE sys- 



tem (2.4) is natural, up to the fohowing: 



a) roto-translations (Euchdean motions) - which can be taken care of by 
setting a base point and an initial moving frame; 

b) constant multiplications of K (similarity transformations) ; 



c) four possible values of a which all satisfy ( |2.3| ) (one for each quadrant: 
meaning that we will have to choose only one, when we perform the 
construction algorithm for isothermic coordinates). 

Proof. Consider the system 



G 



^ VE^/^ ^^^^ 



whi ch, after simplifications, actually becomes equivalent to system (2.4) 



(2.5): 



( g,g.-2gF. + G. -E,- ^E^+2 F^ 
2{G~^)VE 2E.Rf^ 



+K ^sma)H rns(rw) 

' ve ' ' ./^r^ 



^ _ S^ig.-2§F.+G. _ -Ey- §E^. + 2 F. 

' 2iG-^)VE 2E,^^ J 



,1^2 ay - f Oa: . 
+K — cos(a) sinlry 
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Let 



Pi{x,y) 



P2{x,y) 



-Ey — ^ Ex-\-2Fx 



2E 



4^ 



in o r der t o simplify the computation. With these simplifications, the system 
(l2ll)-(IZ5l) becomes 



= ((Pi + Qi) cos(q) - (Pa + Q2) sm(Q)) , 
Kf, = K"" (-(Pi + Qi) sin(a) - (Pa + Q2) cos(a)) . 

or in its equivalent form 



^ = ((A + Oi) cos(a) - (Pa + Qa) sm(a)) , 



= (-(A + Qi) sm(Q) - (Pa + Qa) cos(q)) . 



Since the quantity K is always positive, showing compatibility condition 
K^p = Kp^ is equivalent to showing compatibility condition 



7/3 V"'/;37 

Next, we write the mixed partial derivatives of in the following way 



K 



7/3 



Ml 



sin(Q) cos(a) 



^ p _ FJtp 
Sy-^l E Sa;-^1 a 



G 



^ — COS (q + , sm Of) 



K 



OX ^ 



a p 



9y 



Pa 



E dx-' 



G- 



F2 



sin(Q) cos(a) 



F A. 

E dx 



Pi 



sin (a) 



cos (a) 



+ W^2(a) 



In the above computations, we introduced the quantities 



Wi(a) = Pi(cos(a)),3 - P2(sin(Q)),3 + (Qi cos(a) - (52sin(a))^ , 
W2{a) = — Pi(sin(a))7 — Pa(cos(Q))7 + (— Qi sin(a) — Qa cos(a)) , 
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for simplification. The expression Wi{a) - W2{a) is evaluated below: 

Wi{a) - W2{a) = if I -Pi 



: Sin [a) H sm(a) cos(Q;j 



—K \ P2 \ sin(a) cos(a) + 



F2 



COS (a) 



+K\pA -^cos^(a) + 



'■y E ^ 



E 



F2 



sin(a) cos(a) 



-K P2 



sin(a) cos(a) + 



'■y _B"a: . 2 



4^ 



G 



p2 



sin (a) 



— — sin(Q) cos(a) — ^ sin(a) cos(a) 



(t7*) 



COS (a) + 



-if 



-if 



dy^^/E^ ■ f \ / \ I V^~^ „:„2 



F2 



sin(Q;) cos(a) + 



V 



G- 



p2 



sin (a) 



sin(a) cos(a) 



E 



£2 
E 



sin (a) + ( — . )— p= cos (a) 



- sin (a) 



- sin(Q) cos(a) + siri(a) cos(a) 
E 



-K 



-K 



Tin 

sin(a) cos(a) H 7====^ cos'^ (a) + , 



9h ^ fT, W ' „ 2 



: sin(a) cos(a) 



J 



jr2 



cos (q) + ( — . )— p= sin (a) 



WE 



K 



F2 



a / \ 

arc V Vfi/ , ~ B*^^ \ « 



V 



y/E 



+ ( 



G- 



G- 



F2 



\ 



J 
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Next, we impose the compatibility condition 




J 



After dividing the equation by K and replacing the expressions of Pi and 
P2, we obtain the condition stated in the theorem. □ 



Corollary 2.2. In case F = 0, the necessary and sufficient condition (2.17) 
reduces to 



1 d f Gx 



1 d 



VGdy\2GVEj Ve dx \2^E 



Ey 



GxOlx + EyUy 

2GE 



9 I ax 

9x\^l , aya^ dyy^/G' 



Ve 



+ 



/GE 



+ 



; V .ml 



fG 



Corollary 2.3. For isothermal c oordin ate lines, G = E and F = 0, the 

necessary and sufficient condition (2.17) reduces to 



d_ 
dx 



+ 



- 0. 



Corollary 2.4. In case F = 0, m = and away from umbilic point s, the 
rotation angle a is zero and the necessary and sufficient condition (2.17) 
reduces to 



In 



= 0, 



or in its equivalent form 

G 

for any arbitral smooth functions (j){x) and ip{y). 



Finally, we are resuming the construction of the isothermic coordinates, 
based on the algorithm we presented starting from an initial coordinate chart 
{x,y), followed by a well-chosen rotation of angle a{x,y) and a re-scaling 
K{x,y), all up to roto-translations and similarity transformations. 
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Example 1. How is the actual construction algorithm handled numerically? 
In practice, we integrate along 7 and (3 and draw a virtual mesh to construct 
the surfaces: 

/(7,/3) = /(0,0)+ f f^{Q,P{x,y))dp+ f fp{-i{x,y),P{x,y))d-f. (2.18) 
Jo Jo 

We start constructing the surface in a neighborhood of (0, 0) by solving a 
succession of PDE's. First, we fix j = and move along (3, then we fix 
(3, and move along 7. Because of the compatibility condition, reversing the 
order of integration would not change the value f{^,f3). 



Note that in (2.18) we integrate along (7,/?), while the integrands are 



given in terms of the initial parameters {x,y). In order to integrate, we 
need first to solve and keep track of all the quantities 

x{l,f3), y(7,/3), f-({x{l,P),y{l,P)), and /;3(x(7,/3),y(7,/3)), 



and with these solve for /(7, /3) in (2.18). Thus, the problem is solved in 
two steps: 

1. First, constructing the mapping x(7,/3), y(7,/3) and the vector fields 
f^, fp by solving numerically, with a fourth- order Runge-Kutta method. 



the systems (2. 13) -(2. 16), and (2.1) -(2.2), respectively. 



2. Then, constructing the surface locally by solving numerically the sys- 



tem (2.18) with a fourth- order Runge-Kutta method. 



The Unduloid is an example of a CMC-Delaunay surface which is gen- 
erated by revolving a curve in arclength s, namely C = {X{s),Z{s)) in the 
xz-plane, around the z-axis. This surface of revolution is therefore parame- 
terized as 



f{s,e) 



X{s) cos 9 
X(s)sme 
Z(s) 



where the angle 9 takes values between and 2tt. 

Note that the parameters {s,9) are not isothermal. 

We applied our algorithm starting from the parameterization {s,9) (in 
the role of original chart (x, y) ), and constructed a pair of isothermic coor- 
dinates, (/3, 7). 

The figure below represents the CMC unduloid that we obtained by ap- 
plying our construction algorithm. 
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Figure 1: Unduloid surface in isothermic parametrization. 
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